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Currently, the only known way to obtain experimental information about the space-time structure 
of a heavy-ion colhsion is through 2-particle momentum correlations. Azimuthally sensitive HBT 
interferometry (Hanbury Brown - Twiss intensity interferometry) can complement elliptic flow mea- 
surements by constraining the spatial deformation of the source and its time evolution. Performing 
these measurements on photons allows us to access the fireball evolution at earlier times than with 
hadrons. Using ideal hydrodynamics to model the space-time evolution of the collision fireball, we 
explore theoretically various aspects of 2-photon intensity interferometry with transverse momenta 
up to 2 GeV, in particular the azimuthal angle dependence of the HBT radii in non-central colli- 
sions. We highlight the dual nature of thermal photon emission, in both central and non-central 
collisions, resulting from the superposition of QGP and hadron resonance gas photon production. 
This signature is present in both the thermal photon source function and the HBT radii extracted 
from Gaussian fits of the 2-photon correlation function. 



I. INTRODUCTION 

Ideal hydrodynamic calculations [TJ |5J |3J IH |5J [SJ |7] 
of heavy-ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) generate momentum distributions of par- 
ticles that correctly describe the effects of the collision 
fireball dynamics on the elliptic flow signature and suc- 
cessfully reproduce the main aspects of the measured 
momentum spectra for a large number of hadron species 
[SllgllinillllllllllalllllllS]. The good agreement between 
theory and experiment of the single-particle momentum 
distributions leaves, however, the fireball geometry essen- 
tially unconstrained. To eliminate the geometric ambi- 
guities, spatio-temporal aspects of the reaction must be 
explored |16j . The only known way to obtain experimen- 
tal information on the space-time structure of the parti- 
cle emitting source in heavy-ion collisions is through two- 
particle momentum correlations [T71 [TSl HH] . The method 
of two-particle intensity interferometry, originally devel- 
oped by Hanbury Brown and Twiss (HBT) [20^ to mea- 
sure angular distances of stars and other stellar objects, 
exploits quantum statistical correlations between particle 
intensities in two-particle coincidence measurements to 
access spatial information on the emitting source. Even 
though the space-time picture of the source extracted 
from intensity interferometry is necessarily incomplete 
[TBI [TOl HU [21] , it yields powerful geometric constraints 
which supplement the momentum-space information con- 
tained in the single particle spectra. 

However, most theoretical calculations of the HBT 
radii, expressed through variances (Gaussian widths) of 
the source function, do not correctly reproduce their mea- 
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sured momentum dependence [111 [23] . Various ways to 
explain this "HBT puzzle" and possibilities to correct 
this failure have been suggested. They include a mi- 
croscopic treatment of the final freeze-out stage through 
hadronic rescattering |24| , exploration of fluctuations in 
the initial state [5S], different sets of initial conditions 
leading to stronger longitudinal |26l and transverse |27j 
hydrodynamic acceleration, as well as the inclusion of 
viscous effects [551 HHJ ISHl 131] and final state interac- 
tions [321 1331 1311 135] . The influence on the characteristic 
HBT radii of finite resonance decay lengths after ther- 
mal freeze-out was explored [271 13SI [33 , and the impact 
of non-Gaussian features in the two-particle correlation 
function on the extraction of HBT radii from Gaussian 
fits to the correlator and their discrepancy from radii 
extracted through source function variances were stud- 
ied [3H]- Recent work [31] suggests that a comprehen- 
sive approach that includes all these effects in a theoreti- 
cally coherent fashion may be able to simultaneously de- 
scribe single-hadron spectra and two-hadron correlations, 
lending support to the claim that hydrodynamic models 
successfully describe both the dynamics and space-time 
structure (size, shape, and lifetime) of the collision fire- 
ball at the point of hadron emission. 

An interesting question that cannot be directly ad- 
dressed by measuring hadron distributions is the time- 
evolution of the collision fireball. The reason is that 
hadrons interact strongly and decouple late. The obser- 
vation of strong hadron elliptic fiow [5], i.e. a large mo- 
mentum anisotropy in the final hadron spectra, has been 
linked to early thermalization in the fireball and to almost 
perfect liquid behavior of the quark-gluon plasma (QGP) 
created in RHIC collisions [35] [ID] . This has opened the 
prospect of exploring details of the QGP equation of state 
and of its conversion to hadrons [15] as well as its trans- 
port properties (shear and bulk viscosity) through elliptic 
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flow measurements (HJ |33J |331 Hi] . Since the QGP ex- 
ists only in the early stage of the collision, a more direct 
probe of its properties and dynamical evolution is de- 
sirable rather than through soft hadron spectra emitted 
after its decay. We here discuss the emission of ther- 
mal photons which occurs throughout the fireball evolu- 
tion, but most strongly during its hottest, earliest stage 
^5, 46, 47, 48j 49J. Photons interact only electromagnet- 
ically and thus escape from the collision fireball without 
rescattering. For this advantage over hadrons one has 
to pay with a correspondingly smaller production cross 
section which makes direct photon measurements much 
more difficult than the observation of soft hadrons. 

The possibility of using photon elliptic flow mea- 
surements to access the early evolution of momentum 
anisotropies in the expanding QGP fireball was ex- 
plored in [50l |5T1 [52] . The growth of these momentum 
anisotropies is driven by spatial anisotropies of the pres- 
sure gradients in the initial hot source created in non- 
central heavy-ion collisions. The collective flow response 
of the QGP to these gradient anisotropies depends on 
its speed of sound and its viscosity [53, ^54] . We propose 
that, in a similar way as photon elliptic flow opens a 
window on the early evolution in momentum space, two- 
photon HBT interferometry may help to constrain the 
spatial size and deformed shape of the fireball during this 
stage. By providing access to both the space-time and 
momentum-space structure of the source at early times, 
thermal photons may thus help to constrain the validity 
of hydrodynamic models, the QGP EOS and its transport 
properties during the crucial hottest and densest period. 

While photons are emitted throughout the evolution 
of the fireball and thus do not allow for sharp snapshots 
of the early stage alone, their emission rate depends on a 
high power of the fireball temperature [45", 1^, TT, '45', 15] 
which gives preference to early emission 55j 56, 5Z,i 581 
[59] . Indeed, we will show that the photon emission func- 
tion has a distinctly bimodal structure, with two compo- 
nents that refiect emission of photon pairs of given pair- 
momentum from different spatial regions at early times 
(when radialcollective flow is still weak) and at late times 
(when the flow is strong) .Although a single measurement 
is insufficient to cleanly separate these components, sys- 
tematic exploration of the dependence of photon HBT 
radiion magnitude and direction of the photon pair mo- 
mentum may eventually allow to do so. 

The bimodal structure of the photon source is only 
one of two sources for strong deviations of the emission 
function from a simple Gaussian shape. The other is the 
massless nature of the photon. For two-photon correla- 
tions it is therefore even more important than for hadrons 
that the HBT radius parameters are computed through 
a procedure that matches their experimental extraction. 
We here use a generalization to non-central collisions of 
the Gaussian fit technique developed in (38] [60]. Since 
photons, unlike charged hadrons, do not suffer from final- 
state Coulomb interactions, this technique is even more 
appropriate here than it is for hadron correlations. 



Our study is exploratory in nature and, as such, lacks 
many features needed for a quantitative analysis. To sim- 
ulate the dynamical evolution of the matter created in a 
heavy-ion collision we use AZHYDRO, a (2-t-l)-dimensional 
code developed by P. Kolb [2^ for solving the equations 
of ideal fluid dynamics in two transverse dimensions as- 
suming boost invariance along the longitudinal (beam) 
direction. We use the same photon emission rates as em- 
ployed in earlier studies of thermal photon spectra and 
elliptic flow [501 EU [62] . However, we exclude hard pho- 
tons produced before the flreball medium has thermalized 
in order to focus on thermal processes.. We also assume 
that photons from the post-freeze-out decay of hadronic 
resonances are created sufficiently late that they do not 
correlate with thermal photons and therefore do not af- 
fect the shape of the two-photon correlation function ex- 
cept at immeasurably small relative momenta. (They 
contribute to the single-photon spectra and will thus af- 
fect the normalization of the two-photon correlator, i.e. 
the correlation strength; this increases the statistics re- 
quired for an accurate correlation measurement.) Our 
equation of state assumes chemical equilibrium in the 
hadronic phase below the critical temperature Tc where 
the QGP converts to hadrons, and thus it does not cor- 
rectly reflect its measured |63j non-equilibrium chemi- 
cal composition. All these deflciencies can be removed 
in time before two-photon correlation measurements be- 
come technically feasible. 



II. TWO-PARTICLE CORRELATORS FROM 
NON-CENTRAL COLLISIONS 

A. Correlation and emission functions 

If the emission function is a perfect Gaussian in space- 
time, the two-particle correlation function is a perfect 
Gaussian in relative momentum, and the HBT radii can 
be directly computed from the RMS variances of the 
emission function |22| . Since, however, real emission 
functions are seldom Gaussian, the direct comparison of 
RMS variances with the experimentally extracted HBT 
radii should in general be viewed with suspicion. A more 
reliable (even if more laborious) approach computes from 
the emission function the actual correlation function and 
then performs a Gaussian flt of the latter, using the same 
fit algorithm as employed in experiment |38j . 

The two-particle correlation function for identical par- 
ticles with momenta Pa and pb is defined as a ratio be- 
tween the two-particle coincidence cross section and the 
product of the two single particle spectra as 



C{Pa,Pb) 



p P dN 

TP dN TP 

^^d^^b 



dN 



(1) 



If the particles are emitted independently from the 
source, C{pa,Pb) can be calculated from the single- 
particle Wigner phase space density S{x,K) ("emission 
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function") which describes the probability of emitting 
a particle from space-time point x with momentum K, 
by folding it with the two-particle relative wave func- 
tion [T5] . In the absence of final state interactions (as is 
true for photon pairs) this wave function is a plane wave, 
yielding 



tributions for on-shell particles 



Writing 



Ciq,K) = 
1 

1± 



(2) 
(3) 



1/ d'^xS{x,K) e*9-^| 



9sJd^xS{x,K+l)JdiyS{y,K~ 



where gs is the spin degeneracy of the particles (for pho- 
tons gs = 2), and the -I- (— ) sign applies for bosons 
(fermions). From now on we will only consider the 
bosonic case. 

The correlation function Eq. (|3| depends on the rel- 
ative momentum between the two particle q ~ Pa^Pb, 
qq — Ea—Efj and the pair momentum K = (pa+Pa) 

= {Ea+Eb)/2. For the spectrum Eq. (|2|, the emis- 
sion function S{x, K) must be evaluated on-shell {K ^ 
p), but for the correlation function Eq. ([3| this is not 
necessary [Ull^. 

Since the measured momenta Pa, pt of the particles in 
the correlator are on-shell, the four-momenta q and K 
are necessarily off-shell. For pairs of identical particles 
the relative and pair momenta satisfy the orthogonality 
relation 



(4) 



With this the argument of the plane wave becomes q-x — 
q°t-q-x, with q° = Ea-E,, = fS q, where /3 = K/K° = 
2K/{Ea+Eb). 

Two approximations are often used to further simplify 
Eq. ([3|. The first is the "smoothness approximation" 
which assumes that the emission function varies slowly 
over the momentum range where the correlator deviates 
from unity: 



C{q,K) 



J d'^xS{x,K)e"J-'' 
Jd^xS{x,K) 



(5) 



It is accurate as long as the curvature of the logarithm 
of the single-particle spectrum is small [51]. For thermal 
sources created in heavy-ion collisions this is usually an 
excellent approximation [55] . The second is the "on-shell 
approximation" K° = (Ea+Eb) /2 ^ Ek = Vm"^ + K^. 
It allows us to replace the factor /3 in the relative wave 
function by the pair velocity and to substitute for the 
Wigner densities in Eq. (Isl the classical phase-space dis- 



1 



{Ea + Eb) 




we see that it applies as long as q/{2EK) ^ 1, i.e. as 
long as the source radius is much larger than the Comp- 
ton wave length of the particle pair [22l [65] . In heavy- 
ion collisions this holds for all hadron species (including 
pions) at all pair momenta K, due to their large rest 
masses. For massless photons, on the other hand, the 
on-shell approximation breaks down at small pair mo- 
menta. For m = and K <^ q/2 one obtains instead of 
Eq. ([6| 




Hence, for small- momentum photon pairs (3 — K /K'^ w 
2K/q (which obviously differs from the pair velocity 
which has magnitude c — 1). We will see that this has 
interesting consequences for the structure of the two- 
photon correlation function and the i^T-dependence of 
photon HBT radii. 



B. Gaussian sources and RMS variances 

The main characteristic of the correlation function is 
the g-range over which it decays to 1 . The corresponding 
half-widths of the correlator are called "HBT radius pa- 
rameters" or "HBT radii" . For a Gaussian source they 
can be related to the "homogeneity lengths" (defined be- 
low) of the source [19] . In the absence of final state in- 
teractions they can be extracted by fitting the correlator 
C{q, K) with a Gaussian function of the form 



C{q,K) = 1 + X{K) exp 



E 



q^qjRliK) 



, (9) 



where the correlation strength parameter X{K) is intro- 
duced to account for uncorrelated pairs that contribute to 
the single-particle spectra in the denominator but not to 
the correlated part of the two-particle cross section in the 
numerator. Uncorrelated pairs arise from far-separated 
decay products of long-lived resonances and (if the parti- 
cles carry spin) from pairs with unaligned particle spins 
(e.g. photons with opposite helicity). The correlation 
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strength is also reduced if particles are not emitted inde- 
pendently, but partially coherently |22] . 

In general the HBT radii R^j depend on the pair mo- 
mentum. They are interpreted as width parameters of 
the effective source ("homogeneity regions" [67 ) from 
which particles with momentum K are emitted. In col- 
lectively expanding sources, where the particle momenta 
are correlated with position by a boost with the local 
flow velocity, the homogeneity regions for particles with 
a given momentum constitute only a fraction of the entire 
fireball. However, only the homogeneity lengths can be 
measured with HBT correlations. Furthermore, the HBT 
radii measure only certain combinations of spatial and 
temporal width parameters ( "variances'') of the source, 
as imposed by the mass-shell constraint Q which implies 
q - X = —q ■ {x~f3t): 



((i, -ft?) (5,-/3, t)) 

1 d^C{q,K) 

^2 



dqidqj 



(10) 



Here (/) denotes an average over the the emission func- 
tion. 



_ Sd^xf{x)S{x,K) 
Jd^xS{x,K) 



(11) 



and x^ = x^ — {x^)- The relations ( |10| are exact for 
Gaussian sources where the inverse width of the correla- 
tor agrees with its curvature at q = 0. In this case they 
can be used as a 'shortcut' for the calculation of HBT 



radii, by evaluating the RMS variances in (10 1 directly 



from the emission function S{x, K) instead of first com- 
puting the correlator from ([5| and then fitting it with a 
Gaussian as in ([9]). Here we will not use this shortcut 
but show results obtained with the help of Eq. ( 10 1 only 



for comparison purposes, to illustrate the importance of 
non-Gaussian features in the source and correlators. 



Gaussian fitting procedure for non-central 
collisions 



The emission function S{x, K) is computed from the 
hydrodynamic model AZHYDRD, for hadrons as described 
in [TF and for photons as outlined in [SUl [5T1 [51] which 
folds the thermal photon emission rate with the hydrody- 
namic temperature and fiow evolution. Using this emis- 
sion function in the average (111, we compute the corre- 
lation function Eq. (|5| as 



gs {C{q, K)-l)^ (cos(g ■ x)f + (sin(g ■ x)f (12) 



beam direction {Ki = 0) such that K ■ q = K± qo where 
qo denotes the "outward" component of the relative mo- 
mentum (along the emission direction of the pair in the 
transverse plane). The "sideward" component qs is de- 
fined as the one perpendicular to K± and the beam di- 
rection. 

For collisions between equal-mass nuclei and pairs 
emitted with zero longitudinal momentum in the center- 
of-mass frame, the source is reflection symmetric in the 
longitudinal direction, and the general form (|9| reduces 

to iMiEg 

C{q,K^) = 1 + Ae-(«°«°+«'^=+«'^'+2'°«=«°0 , (14) 

where A, Ro, Rs, Ri and Ros are all dependent on K±. 

For central collisions between spherical nuclei, the 
emission function is azimuthally symmetric and all di- 
rections of K± are equivalent (i.e. the correlator and 
HBT radii only depend on the magnitude of K±). As 
a consequence, the cross-term R^^ vanishes [551 EH]- In 
non-central collisions, the source is deformed and initially 
out-of-plane elongated with respect to the reaction plane 
spanned by the impact parameter b (deflning the x-axis) 
and the beam (defining the z-axis). As a result the cross- 
term is now non-zero, and the 2-particle correlation 



function (12 1 and its HBT radii (14) now depend on the 
azimuthal angle <i> of the emission direction K± . This de- 
pendence is both implicit through azimuthal symmetry 
violations in the emission function S{x,K) with which 
the averages in ( 12 1 are taken and explicit through the 



azimuthal rotation between the reaction-plane coordinate 
system (in which the source has been computed) and the 
osl system defined by the emission direction of the pair 
in which the HBT radii ([ul) are determined [6811691170] : 



qr^ = qo cos $ - sin $ , 
qy = qo sin $ + sin <i> . 



(15) 



The HBT radii are determined by fitting the computed 
correlation function with the functional form (14). To 



where 



this end we generalize the Gaussian fitting algorithm de- 
veloped in j38, i60] to include the Ros as an additional fit 
parameter. (Note that Ros cannot be determined from 
1-dimensional Gaussian fits to slices of the correlation 
function along any one of the osl axes; it requires at 
least a 2-dimensional Gaussian fit in the os plane.) We 
write 

ln[.g,(C(q)-l)] = (16) 
In A - {qlRl + qlRl + qfRf + 2qoqsRls), 

evaluate this expression on a suitable set of q-points q^'^'> 
{k = 1, . . . , N), and define as 



q - x^ (Ea-Eb) t - q^ 



qyV-qiz 



(13) 



with E'^^ = Ej^ ± K ■ q + q'^ /A. We will only consider 
"mid-rapidity pairs" with zero pair momentum along the 



N 

E 

k=l 



In [C(q('=)) - 1] - In A + M{R, q^^)) 



(17) 
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with 



M(fi,gW)=2^)gWi?L+ E (CT ■ (18) 



For the bin error cr^ we use 



C(qW) ™ 1 



(19) 



Minimization of with respect to the fit parameters, 



Sin A 



(i = o,s,Z,os) , (20) 



produces a set of five coupled linear equations that can 
be written in matrix form as 



/3 



(21) 



where a and /3 take the values 0, o, s, Z, os (0 is associated 
with the correlations strength A). The vectors on the 
right hand side have the form 



P — {in X,Rl,Rl,Rf, Rig) , 
_^ln[C(qW)-l] 



V, 



Vos = 



(22) 
(23) 

(24) 



..1 ' 

fc=i '-^fc'' 
w (fe) (fc) 

fc=i ^'^fe'' 



while the symmetric 5x5 matrix T has the components 



T — 



T ■ — — 



^ 1 



E 

fc=i 



«)^ 



(26) 



k=l 



fe=l 



«)' 



with = {o,s,l). As before [38], this set of linear 

equations is easily solved through matrix diagonalization 

of Ta0. 



III. TWO-PION CORRELATIONS IN 
NON-CENTRAL AU+AU COLLISIONS 

As a reference for comparison with the two-photon cor- 
relation functions computed further below, we briefly re- 
view the behavior of two-pion correlations and pion HBT 
radii in non-central Au-|-Au collisions at RHIC energies. 
This analysis complements earlier work ^71j which was 
based on the same hydrodynamical model as used here 



but made use of the shortcut (10 1 to calculate the HBT 



radii directly from the RMS variances of the hydrody- 
namic pion emission function. As explained above and 
numerically studied in [38' , this shortcut becomes doubt- 
ful when the emission function is not well described by a 
Gaussian in space-time. Here, we first compute the two- 
pion correlator (|5| numerically from the hydrodynamic 
emission function STr{x, K) [72] and then obtain the HBT 
radii from a 3-d Gaussian fit to this correlation function 
as described in the preceding Section. For simplicity we 
include only directly emitted pions, i.e. we neglect pions 
from post-freeze-out decays of unstable resonances. 

Figure [Tj shows the azimuthal angle dependence of the 
pion correlation function for peripheral Au-|-Au collisions 
at & = 7 fm, in the limit if _l ^ as we change the direc- 
tion of K±. At K± = the correlation function "sees" 
the entire fireball |71j . Correspondingly, for non-central 
collisions at RHIC energies, the outward radius increases 
and the sideward radius decreases as we move away from 
<& = 0, due to the out-of-plane deformation of the source 
at freeze-out. Beyond $ = 90° these tendencies reverse. 
Through geometric arguments one easily sees that these 
oscillations are seen as a clockwise rotation of the cor- 
relator contours, opposite to the counter-clockwise rota- 
tion of the direction of the transverse pair momentum. 
These same oscillations were observed in previous calcu- 
lations that used RMS variances of the emission function 
as proxies for the HBT radii |7T]. 

Figure [2] compares the azimuthal oscillations of the 
pion HBT radii for 200 A GeV Au+Au collisions at 6 = 
7fm with those of the RMS radii, for two values of K±^ 
(0 and IGeV/c). In addition to the results from the 
3-dimensional Gaussian fit described in the previous Sec- 
tion (labeled "3D" ) we also show, for the purpose of com- 
parison, radius parameters extracted from 1-dimensional 
Gaussian fits to slices of the correlation function along 
the qs and qo axes (labeled "ID"). (Note that the R^^ 
cross term cannot be extracted from these slices, so the 
corresponding ID curves are missing.) In the 3D Gaus- 
sian fits, we restrict the qi fit range to a narrow window 
qi < 9max = 0.02 GeV, in order to avoid distortions in 
Ro arising from strong non-Gaussian tails of the correla- 
tor at larger qi values (see Ref. [3F and Figure [t] below) . 
This is permissible since we found that the shape of the 
correlation function along the qi direction is almost in- 
dependent of 4>, including its non-Gaussian tail. The 
reason is that, even though non-Gaussian effects of the 
source are strongest along the beam direction, they are 
mainly caused by the strong boost-invariant longitudinal 
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FIG. 1; (Color online) Contour plots of the pion correlation function C{qo, Qs, Qi—O) — 1, at $ = 0°, 45°, and 90° in the limit 
A'x 0, for 200 yl GeV Au+Au collisions at 6 = 7fm. When rotating K± counter-clockwise, the correlation function rotates 
clockwise. Successive contours are separated by 0.05, starting from C(0) — 1 = 1 in the center. 
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FIG. 2: (Color onhne) Pion HBT radii Ro, Rs, and R^^ from 
semi-peripheral (6 = 7fm ) Au+Au oscillations at y'sjvjv = 
200 GeV, as functions of the azimuthal emission angle See 
text for discussion of different extraction methods. 



expansion of the fireball which is independent of impact 
parameter and transverse position in the reaction zone, 
and thus ^-independent. If we use a larger fit range in 
qi , thereby capturing more of the non-Gaussian tail, we 
find an intercept parameter A < 1 that oscillates with 
These oscillations interfere with the oscillations of 
the transverse radius parameters Ro and Rs , significantly 
modifying their oscillation amplitudes. This is an unde- 
sirable effect caused by non-Gaussian features which lie 
entirely in the longitudinal domain. By restricting the 
fit range in the qi direction as described we can keep the 
correlation strength parameter close to 1 for all angles 
and ensure that the azimuthal oscillations of the trans- 



verse HBT radii faithfully reflect the oscillations in the 
widths oi C{qo,qs) as seen in the contour plots. 

Figure |2] shows that the azimuthal oscillation ampli- 
tudes of the HBT radii from the Gaussian fit agree very 
well with those of the RMS radii even though (as pre- 
viously observed |38j ) the ^-averaged HBT radii differ 
somewhat from the corresponding RMS radii. R^^ is seen 
to increase with K± , continuing the previously predicted 
[711 [75] and experimentally confirmed [74 trend to larger 
values of Kt- 



IV. PHOTON INTERFEROMETRY 

Since photons are emitted throughout the fireball 
evolution, with emission at higher transverse momenta 
weighted towards earlier times and higher temperatures, 
we expect (and our calculations confirm) that, for suf- 
ficiently high K±, the thermal photon HBT radii will 
reflect geometric and dynamical characteristics of the 
smaller and more deformed early source. This is a simple 
consequence of the dominance of QGP radiation in the 
single photon yield at high p±. However, similar to pion 
interferometry, at any given K± contributions from lower 
temperature regions that are blue shifted by collective 
flow will mix with the early emission contributions. This 
lower temperature emission is generated at later times 
when the source is larger and less deformed. There are 
two competing processes; small radii with weak flow and 
large source deformations reflected in the early emitted 
photons and larger radii with smaller spatial deformation 
which characterizes the emission from the more strongly 
expanding later stages, superimposed by a non-trivial K± 
dependence driven by longitudinal boost dynamics. The 
complexity of these competing processes demands a com- 
prehensive study of two-photon interferometry and the 
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corresponding HBT radii. 

Bass et al. |59j explored photon interferometry for cen- 
tral Au+Au collisions using ideal hydrodynamics coupled 
to a parton cascade to simulate the early preequilibrium 
emission ^75j. They did not, however, investigate the az- 
imuthal behavior of the photon HBT radii in non-central 
collisions. In this work we pioneer the exploration of 
anisotropics of the photon HBT radii from non-central 
heavy-ion collisions, using the same (2-|-l)-dimensional 
hydrodynamic source exploited in our earlier work on 
photon elliptic flow [501 HH ED IS2 to provide us with 
the thermal photon emission function. We examine pho- 
ton HBT radii from both central and non-central col- 
lisions and compare them with the corresponding pion 
HBT radii. 

Before proceeding we must address several important 
differences between pion photon interferometry. First, for 
pions we were able to ignore spin. Photons carry spin 1 
and possess two possible helicity states, and only photons 
with aligned helicities contribute to Bose-Einstein corre- 
lations through wave function symmetrization. Conse- 
quently, the strength of the photon correlator is reduced 
by a factor of 2 [7B] . The correlator between two photons 
with momenta Pa = K+^q and pi, = K—^q, averaged 
over helicities, thus reads 



curvature of the correlation function C{q, K) at the ori- 
gin. Equation (28 1 shows that for sources with non-zero 
emission duration (P) the outward and sideward RMS 
radii (and thus the curvature radii at the origin of the 
correlator in the corresponding q directions) differ from 
each other at all values of K±. This is contrary to the 
case of massive hadrons where the temporal contribu- 
tions to the RMS radii are suppressed at K± — > by 
powers of the pair velocity /? ^ such that generically 
Rs{K±=0) = Ro{K±=0) (for a discussion of possible ex- 
ceptions see [77]). 

The validity of the on-shell approximation is restricted 
to the region q <^ ^Ej^ which for midrapidity photon 
pairs shrinks to zero as K± 0. This implies that, 
for Ki = photons with small transverse pair momenta 
K± < l/i?sourcc (where i?sourcc characterizes the size of 



C(g,K)«l 



J d'^x S{x, K) 



(27) 



the emission region) the RMS radii (10 1 describing the 
curvature of the correlator a.i q — Q can no longer be ex- 
pected to faithfully represent the width of the correlator. 
In other words, for photons the standard connection that, 
for Gaussian sources, relates the inverse width of the 2- 
photon correlator to the size of the photon-emitting re- 
gion is broken for small pair momenta. Consequently, for 
soft photons the HBT radii extracted from a Gaussian fit 
to the 2-photon correlation function cannot be directly 
computed from the RMS widths of the photon emission 
function, even if the latter is perfectly Gaussian. At large 
K± 3> l/i?source this remains possible for Gaussian pho- 
ton emission functions. However, since hydrodynamic 
photon emission functions are not sufficiently Gaussian 
(see below), the RMS shortcut should be avoided alto- 
gether and replaced by a Gaussian fit to the numerically 
computed 2-photon correlation function. 

We illustrate this issue in Fig. [3] by studying the small- 
q behavior of the 2-photon correlation function for soft 
midrapidity photon pairs with K± = 0.01 GeV/c, com- 
puted numerically from the hydrodynamic photon emis- 
sion function. Shown are 1-dimensional slices of C — 1 
along the three axes qs, qo, and qi (from top to bottom), 
setting the two other q components to zero. By plotting 
C{q, K) — 1 logarithmically against qf one easily identi- 
fies non-Gaussian features as deviations from linear be- 
havior. The "Gaussian width parameters" i?f are given 
by the inverse slopes of the Hues. Clearly Rf > R^> B?^ 
in the case shown. One sees that C — 1 is a perfect 
Gaussian in the sideward, but not in the outward and 
longitudinal directions. Deviations from Gaussian be- 
havior along the qi direction should be expected, since 
they are also seen in the 2-pion correlation functions 
[38j and can be traced [Ml [7H] to the boost-invariant 
expansion dynamics of our hydrodynamic source. They 
extend over the entire qi range, i.e. \a[C{q^) — 1] can 
not be approximated by a straight line anywhere (ex- 
cept for qi = 0). This is different for the outward corre- 
lator ln[C(go) — 1] (sohd line in Fig. [s]) which can be 
quite well described by two straight lines with differ- 
ent slopes at low and high q^, with a smooth transition 
radii (lOl and (28l continue to correctly describe the near q^ = 0.0004 (GeV/c)^ = {2K±y^. Apparently, this 



where we again applied the smoothness approximation, 
approximating the source functions at K ±^qiii the de- 
nominator of Eq. ([3| by source functions at K. A second, 
less obvious difference between 2-pion and 2-photon cor- 
relations arises from the fact that photons are massless 
which is addressed next. 



A. Failure of RMS variances for soft photons 

It turns out that for soft photon pairs {\K\ < 
l/-Rsource) both the on-shell and smoothness approxima- 
tions discussed in Sec. |IIA| become problematic. In this 
subsection, we consider these approximations in turn. 



1. On-shell approximation 



Consider the source RMS radii Ro and Rg defined in 



(101 



R„ 



(28) 



where for midrapidity pairs /3 = = K^/K^. With 
the on-shell approximation = (£"„ -I- Ei,)/2 « 
(see Sec. II A) this becomes the pair velocity which for 
massless photons is /3 — 1, independent of transverse 
pair momentum. The on-shell approximation is valid for 
q < 2Kj^, so it holds at q = such that the RMS 
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change of slope is associated with the transition between 
the two limits discussed at the end of Sec. Ill Al Such a 
change of slope is not seen for the 2-pion correlator which 
for our hydrodynamical source is quite well described by 
a single Gaussian in the outward direction . 



FIG. 3: (Color online) One-dimensional slices of the 2- 
photon correlation function C{q,K) for soft photons from 
central 200 A GeV Au+Au collisions, for midrapidity pairs 
with K± = O.OlGeV/c. C — 1 is plotted logarithmically 
against (i = o, s, V) to facilitate visual extraction of the 
Gaussian "width parameters" from the slopes of the curves. 
The blue dotted and brown double-dash-dotted lines are com- 



puted from the RMS variances of the source via Eqs. ( 29 \ and 



Eq. (301, respectively. See text for discussion. 



We can try to understand the behavior of the curves 
in Fig. [3^ analytically by writing (for iiT; = pairs) /3 = 
K^jK^q) = I3{q), with the limits Yi\Rq^2K^P{q) = 1 
and Vmiq^2K^ P{q) — 2K±_/q for massless photons, and 
using this to express the exponent in Eq. ([9]) in terms of 
the RMS variances (10 1, assuming a Gaussian source. In 
the limit of small relative momenta q <^ 2Kj_ we obtain 



f qo = qi = 0, 

J2 q^q,{{3:^-P^{q)i){iJ~PM)i)) = \ ql ((-^o) " 2(i„?) + (P)) , q, = qi = 0, (29) 

which is identical to the RMS variance calculations, setting /? = 1 for photons. For large relative momenta q ^ 2K± 
we find instead 

f qs{S:^s), qo = qi = 0, 

q^q,{{i^-P^{q)i){i,~|3,{q)t)) = I ql{il) - AKMi^oi) + 4^1 (P), q, = qi = 0, (30) 
iJ I qK^^f), qo = qs=0. 



We see that in this second limit factors of i come with- 
out factors of qo, but with factors of K± instead, and 
this changes the g-dependence of the correlator in qo 
direction. We note that, if we had considered photon 
pairs with non-zero pair rapidity, an analogous differ- 
ence would have also appeared between the last lines in 
Eqs. (p9l) for qo < 2K and ([30| for g > 2K, where 
K 



Figure [3] assumes Ki = 0, hence for a Gaussian source 
qgRf — qzl^f) holds at all values of q; the change of slope 
of the longitudinal correlator seen in Fig. |3] thus reflects 
deviations of the emission function from a Gaussian form 
along the longitudinal direction. We will discuss these in 
more detail further below. The sideward HBT radius is 
always free from temporal factors, since Ps = by defi- 



nition. The sideward correlator in Fig. [3] exhibits a per- 
fectly constant slope, indicating that the hydrodynamic 
photon emission function is well described by a Gaussian 
in Xs direction. 

At large relative momentum qo the first term in the 
middle line of Eq. (30) dominates, and the slope of 
ln(C — 1) converges to (Xq), with no dependence on the 
emission duration left. This is interesting because it al- 
lows to separate the geometric from the temporal struc- 
ture of the source with a single accurate measurement 
of the correlator at a fixed value of K± , by studying its 
different slopes at small and large qo- (We are aware, 
of course, of the experimental challenges of measuring 
thermal photon correlations at very small values of the 
pair momentum K±.) With massive hadron pairs this is 
impossible: For them Eq. (28 1 is always a good approx- 
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imation (as long as the source is sufficiently Gaussian), 
so a separation of geometrical and temporal contribu- 
tions requires measurements for different pair velocities 
(3 ~ i^j,, with the additional assumption that the val- 
ues of (i^), {xt) and {t^) don't themselves exhibit strong 
ifj^-dependence. 

According to Eq. ( pO] ), at large qo the outward slice 
of the correlator looks like a Gaussian with slope (i^) 
and reduced intercept A = }^e~'^^^^*' ^. At small relative 



momentum, its slope is characterized by 



(P) > (ig); this is bigger than the large momentum slope 
since for the hydrodynamic source both —2{xot) and (P) 
are positive. Eqs. (29 1 and (30 1 thus explain qualitatively 



the change of slope along the outward direction seen in 
Fig.lD 

The dotted and double-dash-dotted lines in Figure [3] 
explore to what extent they also do so quantitatively. 
We see that they don't: the dotted line computed from 
Eq. (29) for qo ^ 2K is much too steep whereas the 



double-dash-dotted line computed from (30 1 for q ^ 2K, 
while having the correct slope, is much too low. The 
reason for this failure is subtle: In the limit K±^ = 0, 
our photon emission function is independent of space- 
time rapidity rj since we use a hydrodynamic model with 
longitudinal boost invariance. In this limit the variances 
{z'^) and (P) are infinite. For small non-zero K± they are 
finite but very large. When this happens, the expression 
^ , together with Eq. ( 10 ) , can no longer be used because 
it relies on a Taylor expansion of the correlator around 
q = 0, keeping only the first non-trivial term and re- 
exponentiating the result. This works only for Gaussian 
source, but for Kj_ — > our source becomes very non- 
Gaussian along the rj direction, leading to a failure of 
this procedure in the terms containing (z^) and (P). 

The Gaussian approximation continues to work for the 
first term ~ (5^) in the middle Eq. (30) which is why the 
double-dash-dotted line in Fig.[3]correctly reproduces the 
slope of the outward correlator slice (black solid line) . In 
fact, this slope is the same as for the sideward correlator 
slice (red dashed line), indicating (Xq) = (5;^) as should 
be the case for the central collisions studied in Fig. [3] due 
to azimuthal symmetry around the beam axis. 

As a check for the above analysis we performed a calcu- 
lation of the correlator with a modified source where we 
multiplied the photon emission function from AZHYDRD by 
hand with a Gaussian cutoff in rj direction, S{x, K) i— > 
5(x, K) X e~'' . This simulates a fireball with finite rapid- 
ity width A?7 = 1 /\/2. Now the longitudinal and tempo- 
ral variances (z^) and (P) remain finite even for K±^ — s- 0, 
and the non-Gaussian features of the underlying boost in- 
variant hydrodynamic emission function (which set in at 
larger -q values) are suppressed. Fig. [4] shows that in this 
case the slopes and magnitudes of the outward correla- 
tor slice are very well reproduced by Eqs. (29) and (30). 
The problem in Fig. |3] is thus entirely due to the boost 
invariance of our hydrodynamic model. 

As a corollary we note that, in contrast to hadrons 
whose rest mass always restricts their longitudinal ho- 
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FIG. 4: (Color online) Same as Fig. [s] but for a modified 
emission function that has been multiplied by hand with a 
Gaussian cutoff factor e^^ (see text for discussion). 



mogeneity regions to small subvolumes whose longitudi- 
nal extent is controlled by the inverse of the longitudinal 
expansion rate [211 122] , low-momentum photons explore 
the full longitudinal size of the expanding fireball. They 
can thus tell the difference between a longitudinally finite 
and an infinite boost invariant source. Quantitatively re- 
liable predictions of 2-photon correlations for pairs with 
small transverse momentum K_\_ must therefore be based 
on realistic 3-dimensional evolution models that do not 
assume a boost invariant density distribution. 



2. Smoothness approximation 

In addition to the small K± breakdown of the on-shell 
approximation, the smoothness approximation also cre- 
ates deviations which must be examined. The smooth- 
ness approximation is accurate as long as the curva- 
ture logarithm of the single-particle spectra is small [M] . 
However, with the emission rates used here the photon 
spectra grow very rapidly at low momentum (see Fig 1 
of [51,), suggesting that in this region the smoothness 
approximation may also break down. 

To explore the validity of the smoothness approxima- 
tion we write the correlator (|3| as 



C{q,K) = l±-D{q,K) 

9s 



f d^xS{x,K) e' 



Jd'^xS{x,K) 



(31) 



where the last factor invokes the smoothness approxi- 
mation ([5]) while deviations from this approximation are 
captured by the "correction factor" 



D{q,K) 



\jd^xS{x,K)\' 
JdixS{x,K+^)JdixS{x,K-^)' ^ ' 



When the smoothness approximation is valid, _D w 1. 
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relator slices are largest in the Qs direction. They matter 



FIG. 5: (Color online) The smoothness correction factor (321 
along the outward (go, black solid), sideward {qs, red dashed), 
and longitudinal {qi, green dot-dashed) relative momentum 
axes, for central Au+Au collisions at 200 A GeV and for pair 
momentum K± — 0.05 GeV. For comparison the correction 
factor along the sideward direction qs is also shown for pairs 
with larger transverse momentum Ki_ = 0.2 GeV (blue dotted 
line). 



We note that the correction factor D is much harder 
to evaluate than the "smoothed" correlator (|5| which 
can be efficiently computed by Monte Carlo integration 
with S as weight function. For this reason we usually 
evaluated the 3-dimensional correlator C{q, K) in the 
smoothness approximation and used this as the basis for 
our 3-D Gaussian fits of the HBT radii. Only for very 
small K± we also computed the correction factor D (see 
Fig. [s] below). In the present subsection we check the 
validity of the smoothness approximation by computing 
the correction factor D{q, K) only along the qo, Qs and 
qi axes. 

Figure |5] shows the "smoothness correction factor" 
D{q, K) as a function of relative momentum for pho- 
ton pairs with K± = 0.05 GeV. As the relative momen- 
tum increases, the ratio D begins to deviate from unity 
for both the outward (solid black) and sideward (dashed 
red) slices. Along the qi direction (green dash-dotted) 
one sees almost no violation of the smoothness approx- 
imation. The characteristic momentum which sets the 
scale for violations of the smoothness approximation is 
again q = 2K, i.e. the same scale that characterizes the 
breakdown of the on-shell approximation. Due to the 
infinite longitudinal size of our boost invariant hydrody- 
namical source, the single photon spectrum diverges at 
zero photon momentum; as a result, the denominator of 



the ratio D in Eq. ( 32 ) goes to infinity along the outward 
direction when qo = 2K±, generating a zero of D{qo, K±) 
at this point. No such zero occurs in the sideward and 
longitudinal directions where ± ^ never vanishes. 

As we will see below (see Fig. [?]), the effects of the 
"smoothness correction factor" on the 1-dimensional cor- 



only for small K± < 1 / Rs, 



To illustrate this we show 



as blue dotted line in Fig.|5]the correction factor D along 
the qs direction for pairs with K± = 0.2 GeV (which is 
larger than the inverse source radius in central AuH-Au 
collisions) . In this case the deviations from unity are neg- 
ligible over the entire range where C(q, K) deviates from 
unity. 



3. Non-Gaussian behavior along the longitudinal direction 

The non-Gaussian character of the longitudinal slice of 
the two-photon correlator in Fig. [3] deserves some further 
discussion. It appears to be much stronger than previ- 
ously observed for pions: as we will see below, a ID Gaus- 
sian fit to the longitudinal correlation function leads to a 
^ 20% reduction of intercept parameter A/ whereas the 
analogous reduction for the two-pion correlator is only 
about 5% [35]. This suggests that the non-Gaussian ef- 
fects may be to some extent controlled by the rest mass 
of the emitted particles. 

To explore this issue, let us consider the following 
simple model for an expanding thermalized Gaussian 
source |22J: 



S{x,p) ~ if TOj^ cosh(77) exp f — 



X exp 



T 



2Ar 



2(Ar)2 



(33) 



The constant K contains factors such as particle spin, 
fugacity, and the emission duration At. The emission 
function ( 33 1 has been studied extensively (see for a 
review). With our hydrodynamical model it shares the 
exact longitudinal boost invariance. By expanding the 
Boltzmann factor around 77 = 0, 



T 



2 



(34) 



and evaluating the emission function integrals for the 
RMS variances ( 10 1 by saddle point integration one finds 



R 



, T 



(35) 



where tq is the mean freeze-out time. This result suggests 
T jm^ scaling for the longitudinal relative momentum 
(jqi) dependence of the two-particle correlation function. 
To isolate rest mass effects on the correlation function, 
we should study it as a function of the scaled momentum 
qi ^jTjra^. 

Figure |6] shows longitudinal slices of the two-particle 
correlation function for spinless bosons of varying mass 
( "pseudopions" ) , computed from the hydrodynamic 
emission function on the hadron freeze-out surface for 
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FIG. 6: (Color online) Longitudinal slices of the two-particle 
correlation function for spinless bosons ( "pseudopions" ) in the 
limit of zero pair momentum, for different masses as given in 
the legend. The correlation functions were computed from 
the hydrodynamic pion emission function on the hadronic 
Cooper-Frye freeze-out surface, for 200^4 GeV Au-|-Au colli- 
sions at zero impact parameter, by varying the pion rest mass. 
They are plotted logarithmically against the scaled variable 



See text for discussion. 



central Au+Au collisions at RHIC energies. The corre- 
lation functions are plotted logarithmically against the 
scaled variable q? : the non- linear decrease of the cor- 
relators in this representation exposes their non-Gaussian 
character. At large qi the correlation functions are seen 
to drop much more slowly than expected from their (com- 
mon) asymptotic slope at g/ = 0. These deviations from 
a Gaussian shape are strongest for low-mass particles and 
disappear for heavy particles. At infinite hadron mass, 
the correlator approaches a Gaussian function with a 
slope that agrees with the rest mass independent limiting 
slope at = 0. 

Figure [6] explains why we see stronger non-Gaussian 
effects in the longitudinal correlator for massless photons 
than for pions. Spot checks have shown that the non- 
Gaussian aspects are further enhanced in volume emis- 
sion (as is the case for photons) relative to Cooper-Frye 
surface emission (as simulated in Fig. [6] for "pseudopi- 
ons" ) . Since computing the volume emission function for 
photons is numerically more demanding than the surface 
emission function for pions, we do not present a system- 
atic study of rest mass effects on volume emission. 



B. Photon HBT radii 

In the remainder of this work we will concentrate on 
the transverse HBT radius parameters for thermal pho- 
tons. But since the 3-dimensional correlation function is 
non-Gaussian along the longitudinal direction, this will in 
general contaminate the transverse HBT radii extracted 



from a 3-d Gaussian fit even if the correlator is Gaus- 
sian in the transverse directions [38]. We deal with this 
problem as follows: 

Figure [7] shows 1-dimensional slices of the 3- 
dimensional photon correlation function C(q) — 1 along 
the qo, qs, and qi directions, together with 1-dimcnsional 
Gaussian fits using a ID analogue of the Gaussian fitting 
algorithm described in Sec. II C The fits (solid lines) were 
done in the range < qi < 0.2GeV/c. One notes the 
bad quality of this fit for the longitudinal shce (Fig. [tJI.), 
arising from the strongly non-Gaussian shape of the cor- 
relator in this direction, and the associated significantly 
20%) reduced intercept parameter. First principles 
tell us (and Fig. |6] confirms) that the correlator must be 
Gaussian near qi = 0. The dashed line in Fig. [7] shows 
that by restricting the qi fit range to qi < 0.02 GeV/c (i.e. 
a 10 times smaller window) we can largely cut out the 
non-Gaussian effects and obtain a fit with a reasonable 
intercept parameter close to 0.5 which is not in conflict 
with the similar intercept parameters preferred by the 
ID Gaussian fits in qg and q^ directions. 

The dot-dashed lines show slices of the thermal photon 
correlation function Eq. ^ which is obtained by multi- 
plying the numerically computed curves in Fig.|7^, which 
use the smoothness approximation Eq. ([5|, by the cor- 
rection factor Eq. 
sideward slices are 
"smoothed' 
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Although both the outward and 
lighly sensitive to deviations from the 
correlator at low relative momentum, the 
characteristic momentum q = 2K is sufficiently in the 
tail of the correlator that only the sideward slice exhibits 
these deviations. The deviations along q^ only drive the 
correlator to zero faster and only in the tail when the 
correlator is sufficiently small. 

In the following, we will therefore show transverse HBT 
radii that have been obtained from 3D Gaussian fits to 
the numerically calculated photon correlation functions 
C{q) using a restricted qi fit range \qi\ < 0.02 GeV in the 
longitudinal direction but including in the transverse di- 
rections the full range of qo, qs values where C deviates 
significantly from unity. This procedure is necessary if 
one wants to obtain meaningful values for the azimuthal 
oscillation amplitudes of the transverse HBT radii in non- 
central collisions (see Sec. IVD2I. If these are extracted 
from 3D fits with unrestricted qi range, the non-Gaussian 
features in qi lead to azimuthally oscillating intercept pa- 
rameters A(<I>) which contaminate and unphysically dis- 
tort the azimuthal oscillations of the HBT radii in ways 
that do not reflect the azimuthal dependence of the size 
and shape of the emission function. 

From the 3D Gaussian fits with restricted qi range we 
also extract longitudinal HBT radii. Obviously, these 
reflect the curvature of the longitudinal correlator near 
qi = and do not necessarily provide an accurate de- 
scription of its longitudinal width. The latter is more 
accurately given by ID Gaussian fits along g/, and we 
will therefore show these ID Ri values for comparison. 
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FIG. 7: (Color online) Slices of the photon correlation function C{q) — 1 for central Au+Au collisions at ^/sNN ~ 200 GeV, 
for photon pairs with K± = 0.05 GeV. Figure (a) shows slices of the numerically calculated 3-dimensional hydrodynamic 
correlation function along the qo (circles), qs (squares) and qi (triangles) directions. Note the strongly non-Gaussian shape of 
the qi slice. The other three panels (b-d) show the individual slices superimposed with 1-dimensional Gaussian fits, using a 
fit range qi < Qmax with Qmax = 0.2GeV/c (solid lines) and (in the longitudinal direction only) Qmax ~ 0.02GeV/c (dashed 
line). The dot-dashed line is the full correlation function taking into account the correction factor D in Eq. (32 1 (see text for 
discussion) . 



C. Central collisions 

Figure [SJd shows photon HBT radii from the hydro- 
dynamic model for central 200 A GeV Au-|-Au collisions 
as a function of photon pair momentum K^, obtained 
from Gaussian fits to the numerically computed correla- 
tion function as described above. In out- and sideward 
directions the ID and 3D HBT radii are almost identical, 
except at small K±, reflecting an approximately Gaus- 
sian shape of the correlator in transverse directions. As 
discussed in Sec. IV A| the outward two-photon correla- 
tor develops non-Gaussian features at small K±, and the 
strong non-Gaussian features in the longitudinal direc- 
tion (visible in the large discrepancy between the ID and 
3D longitudinal fit radii Ri) get also worse at small K±. 
Both effects drive the common intercept parameter A for 
the 3D fit (squares in Fig. [sja)) below 0.5, especially at 
low K± , in spite of the restriction of the longitudinal fit 
range to qi < qf^'^^ = 0.02 GeV. For the ID longitudinal 
fit which uses a 10 times larger longitudinal fit range, A 
shows a nearly 25% deviation from the expected value 
of 0.5 as K± 0. This is our motivation for restrict- 
ing the longitudinal fit range in the 3D fits, in order to 
minimize contamination of the transverse radii extracted 
from the 3D fit that stems from a reduction of the com- 
mon intercept parameter exclusively driven by strongly 
non-Gaussian features along the longitudinal direction. 

The dotted curves in Figures [sjjc) (ID radii) andjsjd) 
(3D radii) show the changes to the Gaussian fit radii 
when the "smoothness correction factor" ^ is applied 
to the correlator, as shown in Eq. (31 1. For the ID fit 



effects for all three extracted radii at small K± where 
the factor D in Eq. ( 32 1 introduces non-Gaussian de- 



radii the smoothness correction factor has almost no ef- 
fect along the outward and longitudinal directions; in 
the sideward direction, the broadening of the correlator 
by the smoothness correction factor, shown as the dash- 
dotted line in the lower left panel of Fig. [t] leads to a 25% 
decrease of the ID sideward radius Rs at K± = 0.05 GeV. 
The 3D fits, however, show large smoothness correction 



formations also in the side- and outward directions (see 
Fig. |5]). At Kj_ = 0.05 GeV these, together with the in- 
herent non-Gaussian structure in (/(-direction already at 
the smoothness-approximated level, cause a reduction of 
the common 3D intercept parameter by almost 20% be- 
low its ideal value of |; this, in turn, further distorts the 
three Gaussian fit radii. These distortions happen only 
at very small pair momentum; already at Kj_ =0.1 GeV 
the smoothness-correction effects are almost completely 
gone. 

The most striking feature of Fig. ^h) is the strong 
divergence between Rg and Ro as K± 0. This is due 
to the temporal contributions (in particular the emission 
duration contribution) to i?~ which, for massless photons, 
are not suppressed by pair velocity factors in the limit 
K± ^ as is the case for hadrons. As a result, for 
soft photon pairs the outward HBT radius exceeds the 
sideward one by more than 50%. We have checked that 
this effect disappears, i.e. lim/f^^o(^o — Rs) = for the 
3D Gaussian radii, if we give the photons a small nonzero 
mass. 

In Figure |9] we separate the photon emission function 
into two contributions, corresponding to radiation from 
the quark-gluon plasma and hadron gas phases. Not sur- 
prisingly, the HG HBT radii indicate a larger source for 
the hadronic phase than for the QGP phase whose size is 
reflected in the QGP HBT radii. The measured photon 
HBT radii (QGP+HG) compromise between these val- 
ues, reminding us that photons are emitted throughout 
the fireball evolution. At large K±, the transverse HBT 
radii from the total source approach those from the QGP 
phase, reflecting the fact that at large K± QGP radiation 
dominates the single photon yield. The transverse pho- 
ton HBT radii from the HG phase are not too different in 
magnitude from the corresponding pion radii, except for 
the emission duration effects in Ro at small Kj_ for pho- 



13 




— i?„ ID 


• ID 


ID 


— i?<,3D 


•— • 3D 


R, 3D 



0.0 



0.2 



0.4 0.6 
if^ (GeV) 



0.8 



1.0 



(c) 


— w/o Z? 


\N/D 


(d) 



















• ■ . 






• 





0.05 0.10 0.15 0.200.05 0.10 0.15 0.20 
(GeV) (GeV) 




0.5 1 1.5 
(GeV) 



— HG+QGP 

G«)QGP 
MUG 

200 A GeV Au+Au 
b=Ofm 



FIG. 9: (Color online) 3D Gaussian HBT radii as a func- 
tion of transverse pair momentum for the photon emission 
functions describing radiation from the quark-gluon plasma 
(QGP), from the hadron gas stage (HG), and the total emis- 
sion. 



ent at early time [32 E7j when most of the QGP radia- 
tion is emitted (see also Eq. ([35|), and acerbated by the 
strong non-Gaussian nature along the longitudinal direc- 
tion. (We note that for photons we find Ri < Ro at all 
values of K± whereas the inverse relation Ro < Ri holds 
for pions [38].) 



FIG. 8: (Color online) A"x -dependence of correlation strength 
parameter A (a) and HBT radii (b-d) from 1- and 3- 
dimensional Gaussian fits of the hydrodynamic 2-photon cor- 
relation function for central Au-|-Au collision at RHIC en- 
ergies. For the 3D fits we used g^J" = 0.2 GeV and 
ginax ^ 0.02 GeV; the ID longitudinal' fits were done with 
ginax ^ 0.2 GeV. (b): HBT radii from ID and 3D Gaussian 
fits to the smoothness-approximated correlator ([5|. (c): Low- 
K± blowup of the ID fit radii from (b) (dashed lines), com- 
pared with ID fit radii for the smoothness-corrected correla- 
tor (311 (dotted lines), (d): Low-i^x blowup of the 3D fit 
radii from (b) (solid lines), compared with 3D fit radii for the 
smoothness-corrected correlator (31 1 (dotted lines). 



tons. The difference between Ro and Rg at small K± is 
smaller for the QGP than for the HG emission function, 
indicating shorter effective emission durations for QGP 
photons compared to hadronic photons. 

The longitudinal photon HBT radius from the HG 
phase is significantly smaller than the pionic one, but 
this is at least partially due to the stronger non-Gaussian 
effects in the longitudinal correlation function for pho- 
tons and our restriction of the qi fit range in the 3D fits. 
The QGP photon radii are significantly smaller than the 
HG radii, especially for the longitudinal radius [55], re- 
flecting the strong longitudinal expansion velocity gradi- 



D. Non-central collisions 

1. Photon emission functions 

The results from the previous subsection for central 
collisions have shown that the photon HBT radii can only 
be properly understood if one keeps in mind that photons 
are emitted from all stages of the expanding fireball. This 
is even more true for non-central collisions and for the 
azimuthal oscillations of the HBT radii in this case. In 
this subsection we show that at nonzero K± ~ 1 GeV, the 
photon emission functions exhibit qualitatively different 
characteristics from those we are familiar with from pions 

Figure [TO] shows contour plots of the normalized emis- 
sion function S{x,y; K±), integrated over z and t, for 
midrapidity photons at K± = 0.05 and 2GeV/c. Since 
K± ~ photons are emitted from essentially everywhere 
in the source at any time, but the space-time volume 
covered by the hadronic phase exceeds that of the QGP 
at RHIC energies [5(1, the contour lines in Figure 10 1 
indicate the overall size and shape of the momentum- 
integrated photon emission function, with a stronger 
weight on the late evolution stages. The slight distortion 
of the contours towards the right (positive x direction) 
reflects the nonzero pair momentum K± = 0.05GeV/c 
in <i> = direction. Emission of such photons is domi- 
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FIG. 10: (Color online) Photon emission function S{x, y) 
at A'x = 0.05 GeV « 0, "1> = (upper left panel) and at 
K±^ = 2.0 GeV for three emission directions, "1> = 0°, 45°, and 
90° (other three panels), for 200 A GeV Au+Au collisions at 
6 = 7 fm. The emission functions are normalized to 1 at their 
maxima; contours are drawn in 10% intervals from this max- 
imum value. 



nated by cells which flow hydrodynamically in -fa; direc- 
tion with a flow velocity that has just the right magni- 
tude to boost the peak of the Boltzmann distribution to 
this nonzero K^^ value. Due to the dominance of the HG 
phase for the emission of low-ii'j_ photons [SHI IS] , this 
reflects the strong radial flow in the late hadronic stage 
of the fireball. 

Photons with larger K± — 2 GeV/c, on the other hand, 
are mostly emitted from the early QGP stage [501 ED- 
At this point, radial flow is weak, so we should naively 
expect the photon emission function to peak near the 
center of the flreball, with little distortion by radial flow. 
This is in strong contrast to pions which are only emitted 
at the end of the evolution from the hadronic freeze-out 
surface and thus always feel the full radial flow of the late 
hadronic fireball. For large K± values, hadronic emission 
functions are strongly peaked near the edge of the fireball, 
forming narrow crescent-shaped slivers that straddle the 
surface of the fireball [7T] since it is in these regions that 
one finds the strongest radial fiow which most efficiently 
boosts the local thermal distribution with temperature 
Tf towards large K± values. 

The three Kjl = 2GeV/c panels (Fig. [lO];b-d)) show 
something quite different: For large K±, the photon 
emission function peaks near the center of the fireball (at 
early times, as studies of its t-dependence have revealed), 
but it exhibits a "bulge" in the direction of the emit- 



ted photons that reflects a contribution from later times 
when the matter is in the hadronic phase and boosted 
by strong radial flow. Even though the hadronic phase 
radiates at much lower temperature ((Thg) ~ 150 MeV 
while (Tqgp) is almost twice as high), which strongly 
suppresses its photon emission rate, the larger space-time 
volume covered by the HG phase partially compensates 
for that loss in rate and makes hadronic photon emis- 
sion sufficiently competitive at K± = 2GeV/c to give a 
visible ffow-boosted contribution to the photon emission 
function. The hadronic component is not strong enough 
to collectively "squeeze" the photon emission function 
towards the edge of the fireball (as was the case for pi- 
ons), but it distorts it visibly into the direction of the 
emitted photon pair. This parallels the study of photon 
elliptic flow which showed [50j .51J that, even though at 
K± ~ 2GeV/c the single-photon yield is dominated by 
QGP radiation, the stronger elliptic flow of the HG pho- 
tons still controls the total photon V2, and one must go 
to significantly larger transverse momenta before the el- 
liptic flow of photons provides an uncontaminated view 
of the earliest QGP stage. 

Were it not for the hadronic "bulge" , the K± = 



2GeV/c emission functions shown in Fig. 10 would in- 
deed provide a clean reflection of the original source ec- 
centricity in non-central collisions. The emission function 
in the Fig. 10 >d (dominated by QGP radiation) is clearly 



more eccentric than the one in the Fig. 10 1 (which is dom- 



inated by hadronic radiation). What Fig.jlOjteaches us is 
that the full story is more complex than suggested by this 
naive view, and that any interpretation of azimuthal os- 
cillations of the photon HBT radii must properly account 
for the two-component structure of the photon emission 
function which is generated by the interplay between non- 
boosted QGP and flow-boosted hadronic photon emis- 
sion, with relative weights that depend on K±. 



2. Azimuthal oscillations of photon HBT radii 

Figure [TT] shows the azimuthal oscillations of the 
transverse photon HBT radius parameters in semi- 
peripheral Au-I-Au collisions at RHIC. The HBT radii 
were extracted from a 3D Gaussian fit to the corre- 
lation function, with restricted longitudinal fit range 
Qo < 0.02GeV/c as explained above. Qualitatively, the 
features seen in Fig. 11 agree with those for pions shown 
in Fig. [2j On a more quantitative level, one observes 
important differences: the sideward radius is smaller for 
photons, reffecting preferred emission from the smaller, 
hotter fireball interior and earlier times, while the rela- 
tive oscillation amplitude is larger, refiecting the stronger 
source eccentricity at early times. The ifj^-dependence 
of the azimuthally averaged Rg value is weaker for pho- 
tons than for pions; the likely reason is weaker radial flow 
at earlier times when a large fraction of the photons are 
emitted. The /Cj^-dependence of Rg is much stronger for 
photons than for pions, mostly due to the much larger 
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FIG. 11: (Color online) Azimuthal oscillations of the 
transverse photon HBT radius parameters Ro (upper left), 
Rs (upper right), and R^s (lower left), for photon pairs 
with K± = 0.2,0.6, and l.OGeV/c from semi-peripheral 
200 ^GeV Au+Au collisions at fe = 7fm. 



FIG. 12: (Color online) Normalized oscillation amplitudes of 
the photon (left column) and pion (right column) HBT radii 
from 200 A GeV Au+Au collisions at 6 = 7 fm, as a function 
of transverse pair momentum. 



photon Ro radii at small K± resulting from the emission 
duration contribution which at low K± is suppressed for 
pions but not for photons (see discussion in Sec. IV A I. 

Interestingly, the R^^ cross term for photons oscillates 
with an amplitude that is almost independent of K± and 
even bigger at small K± than at large K±. Fig. [TT] : 
exhibits almost pure sin(2<i>) oscillations for the R^^. The 



relationship [551 171] 

Rl^ = cos(2$)(iy) +sin(2$) 



if 



(36) 



/3(g) ((£t) si 



sin (f> - (yi) cos <i> 



which holds for emission functions whose dependence on 
x,y,t is Gaussian and should thus be reasonably accu- 
rate for our source, indicates that the larger sin(2<i>) os- 
cillation amplitude for photons at low Kj_ is caused by 
a larger source eccentricity for low-momentum photons 
than for low-momentum pions. This is due to the signif- 
icant weight of early photon emission even at zero trans- 
verse pair momentum. 

For pions it was found that Kj_ = pion pairs are 
emitted from the entire fireball, and that the normalized 
azimuthal oscillation amplitude of the sideward radius 
Rs can thus be used to measure the spatial eccentric- 
ity of the fireball at pion freeze-out \7W\. In the same 
spirit we can try to extract the fireball eccentricity from 
the normalized azimuthal oscillations of the photon side- 
ward radius at Kj_ . The sideward radius is the signature 
of choice since it is a purely geometric observable, un- 
contaminated by temporal contributions (see Eqs. (29) 



and (30l). The discussion above indicates, however, that 
photons even at K± — are emitted more often from the 
hot fireball center at early times than from the cooler 
periphery at later times so the normalized azimuthal Rg 



oscillation for photons will measure the effective fireball 
eccentricity at earlier times than for pions. This is, of 
course, exactly what we hoped to obtain. The only un- 
expected aspect is that this would work even for K± — 
where we did not anticipate early emission to play quite 
as important a role as we now see. 

The contour plots in Fig. [TO] show that the space-time 
character of photon emission differs from pions even more 
strongly at larger K± where photon emission is even more 
strongly concentrated at early times close to the fireball 
center whereas pion emission is almost completely sur- 
face dominated and concentrated to a thin sliver near the 
fireball edge f7T| . This makes a geometric interpretation 
of the normalized R^ oscillation amplitude in terms of 
spatial eccentricity of the photon emission function pos- 
sible even at larger K± values. By going to larger K±, 
we can thus measure with photons the fireball eccentric- 
ity at even earlier times, by analyzing their normalized 
Rs oscillation amplitudes as a function of K±. (For pi- 
ons the straightforward geometric interpretation of this 
observable is lost for > [79|.) 

In Fig. [12] we show the normalized oscillation ampli- 
tudes [70 which, following Retiere and Lisa [TS], can be 
expressed through the following normalized 2""^ order az- 
imuthal Fourier components of the HBT radius parame- 
ters i?? ($) [50]: 



n2 



?o,.)(0)-i??o,.)(f) 



Rm 



f i?2(f) 



7?2 



i?2(0) + i?2(|) 



(37) 



The oscillation amplitudes for R^ (Fig. 12 1) and R^ 
(Fig. 12 d) differ mostly by an overall sign for both pho- 
tons and pions; any additional differences in magnitude 
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are due to azimuthal oscillations of the temporal contri- 
butions to i?o [Ml E] ■ The latter are larger for pions than 
for photons. Fig. [T2| ; shows that the azimuthal oscilla- 
tions of the R^^ cross-term are much larger for photons 
than for pions (as already seen in Figs. 11 and[2|, and 
that this difference is mostly due to the earlier emission of 
the photons (it is more pronounced for the photons from 
the QGP phase than from the hadron phase, but even in 
the HG phase this cross-term oscillates more strongly for 
photons than for pions at large Kj_ where HG photons 
are emitted earlier than pions). 

The most important information can be extracted from 
the plo t utilizes the oscillation amplitudes of . Fig- 
ure [l2p shows that R^ oscillates more strongly for pho- 
tons emitted from the QGP stage than from the hadron 
phase. This reflects the larger eccentricity of the fire- 
ball at earlier times. The oscillation amplitude for the 
total photon emission function interpolates between the 
QGP and HG limits: for K± > 2GeV/c it coincides with 
the QGP curve (indicating complete QGP dominance for 
Kj_ > 2GeV/c), but even at ifj, = the QGP contri- 
bution still plays a significant role. Retiere and Lisa [75] 
showed for pions that one can extract the source eccen- 
tricity at freeze-out from the relation (see also Eq. (2) in 



2 

s,0 



(38) 



in the limit K± ^ 0. In Ref. [74] this procedure was suc- 
cessfully apphed to RHIC data. We have argued above 
that for photons a geometric interpretation of the r.h.s. 
in ( 38 1 should remain possible at non-zero K± . If this 



is true, we can follow the fireball eccentricity backward 
in time by following the black solid line in the upper 
left panel of Fig. [12] from low to high K±: K± — 
would represent a time somewhere around whereas for 
K± > 2GeV/c we would probe times close to thermal- 
ization of the QGP. Of course, fixed K± values cannot be 
mapped one-to-one to sharply defined emission times - 
the intention of our argument is to point out a qualitative 
and novel tendency that becomes accessible with photon 
interferometry. The analogous procedure for pion HBT 
oscillations has no similar meaningful geometric interpre- 
tation. 



V. CONCLUSIONS 

In this paper we presented a first comprehensive study 
of two-photon correlations and the azimuthal oscillations 
of the photon HBT radii in non-central heavy-ion col- 
lisions. We based our investigation on the hydrody- 
namic photon emission function for Au-|-Au collisions at 
^/s = 200 A GeV with impact parameter 6 = 7 fm as a 



model system. By comparing photon and pion HBT cor- 
relations we were able to identify several key differences, 
both in the formalism and in the numerical results, that 
make azimuthally sensitive photon interferometry an ex- 
citing prospect for future experimental studies. 

The two most important insights from this study are (i) 
that since real photons are massless, a qualitative change 
of the shape of the 2-particle correlator C{q) for soft pho- 
ton pairs with K < l/i?firebaii develops which, if it can 
be measured, will allow to separate the spatial and tem- 
poral aspects of photon emission in a model-independent 
fashion, and (ii) that at all transverse momenta within 
the range of validity of the hydrodynamic model photon 
emission is essentially volume dominated, with a strong 
QGP component from early times close to the fireball 
center and rather weak distortions from a secondary fiow- 
boosted hadronic emission component that extends more 
towards the edge of the fireball in the emission direction. 
The first observation requires measuring 2-photon corre- 
lations at very small pair momentum which will be very 
difficult; it also requires a careful shape analysis of the 
2-photon correlator in relative momentum since both the 
on-shell and smoothness approximations break down in 
this K± range, and the correlator becomes strongly non- 
Gaussian. The second feature is more useful in practice 
since it works over a large range of pair momenta; it al- 
lows to use photon HBT radii for measuring the size and 
shape of the fireball at early times, before hadronization 
and hadronic freeze-out, and to map out the time evolu- 
tion of its spatial eccentricity, by studying the normalized 
azimuthal oscillation amplitude of the sideward radius as 
a function of transverse photon pair momentum. This re- 
alizes a long-held dream for using 2-photon correlations 
as a microscope for measuring the early fireball geometry. 
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